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We examine the interplay of surface and volume effects in systems undergoing heat flow. In 
particular, we compute the thermal conductivity in the FPU jS model as a function of temperature 
and lattice size, and scaling arguments are used to provide analytic guidance. From this we show that 
boundary temperature jumps can be quantitatively understood, and that they play an important 
role in determining the dynamics of the system, relating soliton dynamics, kinetic theory and Fourier 
transport. 



Transport in low dimensional systems has enjoyed renewed interest due to new efforts to understand why some 
theories have a bulk limit for transport coefficients and others do not The FPU /3 model, in particular, provides 

a classic example of a system with no bulk limit |4[|^. Steady state transport is typically generated by imposing 
boundary conditions on surfaces of the system of interest, that provide the necessary behavior, such as shearing, 
heat flow and so forth. A physical and ubiquitous feature of these boundary-driven non-equilibrium steady states 
are surface discontinuities in certain observables In fluids that are being sheared, a slip velocity develops at 
the boundary reflecting the fact that the fluid just inside the system does not move with the same velocity as the 
walls. In systems undergoing heat flow, as in the FPU model we study here, a 'surface resistance' or boundary jump 
in temperature develops. We will show that these surface effects can be quantitatively understood, and that the 
interplay of surface and volume effects can be used to define dynamical regimes, characterized by different scaling 
properties. To provide a more complete picture, we compute both the temperature and size dependence of the thermal 
conductivity k, as well as the speed of sound and other kinetic theory quantities. From the behavior of k, we will 
also see that the presence of a constant thermal gradient does not guarantee that Fourier's law, specifically the linear 
relation between the heat flow J and VT, J = — kVT, is valid. (We refer to the relation J = — kVT as Fourier's law, 
even in the absence of the bulk limit.) The role of solitons can also be understood. 

Heat flow is generated by applying hot (Tj) and cold (T") sources to the different boundaries of the system. 
However, in general there will be a mismatch between the surface temperature and temperature just inside the 
system, T^, denoted 5Ti = \Ti — T^\ Consequently, the input temperature difference, AT — — T°, can be 
decomposed into the two surface jumps, STi, ST2, and the gradient contribution: AT — J dxVT + 5Ti + 5T2. The 
jumps 5Ti can be expressed as 5Ti = rfUT\i, where VT^ is the internal temperature gradient extrapolated to the 
surface z, and 77 is on the order of the mean free path at that temperature \i 0,^. Hence we write STi = a\i VT|j, 
where a is a constant corresponding to the efficacy of the boundary thermostats as they couple to the system. In 
this work, it will turn out that the boundary jumps become important in the near equilibrium regime, so that they 
are symmetric, 5Ti ~ 5T2, and the gradient is constant, VT^ = VT. Thus we are led to the general near-equilibrium 
relationship, using the length L between the thermostats, 

[2a\ + L] VT AT. (1) 



Seemingly simple, this relation can be used to derive many properties of the FPU model, and others in general. 



In this letter we examine the FPU [3 Hamiltonian, H = ^ + \Tnuj'^{qk+i ~ ^fc)^ + f (<?fc+i — QkY 

the rescaling pk = pJ^w^Vm^, = q'j.LO\/m, we obtain the conventional form of the FPU [i model, 



k=l 



Under 



(2) 



where H/^ — H / (m'^uj^). It is important to note that the temperature and the coupling are noi independent parameters, 
but changing the temperature is equivalent to changing the coupling /3. This can be seen by the further rescaling 
p'k — Pk/\/P, q'k = Ikl^fP, which leads to a unique, dimensionless, Hamiltonian H = Hp^i = ^Hp. Since we have 
rescaled the momenta p\ — Pp'^ , the temperatures in the two formulations H and Hp are related by T = (iT' . So 
the most general physics can be obtained by using Hp around a fixed temperature (for instance T' — \) but with 
different couplings /?, or by using H at various temperatures T . In this study we use H and probe temperatures from 
10"'* to 10'*, where T^ ~ {pD (ideal gas thermometer), using L sites and either free {qo = qi, qL = qL+i) or fixed 
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{qi — Ql — boundary conditions. Perhaps surprisingly, these conditions have an effect on the measured transport 
properties. 

Heat flow is generated when thermostats are coupled to the boundaries of the system. The system evolves according 
to Hamilton's equations of motion except at the boundaries, where we use Nose-Hoover (NH) thermostats or demons 
I pot (see Q, Eqs. (3)). The NH thermostats are observed not to thermalize the boundaries for large temperature 
gradients, a problem particularly acute for small lattices. The demons, on the other hand, do not suffer from this 
problem and hence provide a resolution to the problem recently studied in [ pT[ . We have considered a range of 
coupling strengths as well as thermostatting one or more sites on each edge. Irrespective of the form of the boundary 
thermostats, we note that their introduction breaks global momentum conservation, Ptot ^ 0, while it is conserved on 
average; (Ptot) = 0. The equations of motion, on the other hand, are still invariant with respect to the translations 
Qk ~^ Qk + c for free or periodic boundary conditions, where c is an arbitrary constant. This disparity arises from the 
fact the the equations of motion for the whole system are no longer Hamiltonian. 
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FIG. 1. Size dependence of the thermal conductivity of the FPU model for free or fixed boundary conditions with different 
thermostats. The Nose-Hoover thermostats do not thermalize the boundaries for small L and hence their results deviate from 
those of the demons. L" ''^ behavior will be discussed along with Fig. 0. 



The definition of heat current J depends on the form of the discrete gradient in the Hamiltonian and continuity 
equations. For instance, we may choose hk — -I- U{qk+i — qk) (where H — /i^, U{x) = /2 -\- and 
hk + Jk+i — Jk = 0, leading to the consistent form Jk = —PkU'{qk — qk-i)- Regardless of the form chosen, the 
measurement of J is the same. We then use Fourier's law, J = —k{T, L)VT, to obtain k. In Fig. ^ (NH-fixed), we 
reproduce the analysis of plotting k using the boundary temperatures (r{',T2) = (2.4,15.2) to obtain k(8.8,L) 
and compare this to the results obtained using demons and also the free boundary condition results using the same 
(T°,r2). This figure raises several significant points. First, the NH results do not thermalize the boundaries for 
smaller L, which accounts for the disagreement with the demon results. There is also a clear difference between 
the fixed or free boundaries. Free boundary conditions are found to result in the diffusion of the entire system in a 
Brownian manner: (g^) ~ t, an essential difference. The differences in k disappear as we increase the system size since 
the temperature gradients become small and we reach the near equilibrium regime. For small L, while VT ^constant, 
the system does not obey Fourier's law. We return to this in Fig. || below. Consequently, referring to these points as 
K is not entirely correct. 
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FIG. 2. (Top) T and L dependence of k. One can see two regimes. For large T, there is an approximate T^/'^ behavior 
which can be understood from scaling arguments. For low T, a 1/T dependence which arises due to the weak coupling limit. 
The region where T ~ 10 has k largely independent of T, so that Fourier's law can be strongly violated while the temperature 
profiles remain linear. This region was studied in [5] and Fig. 1. (Bottom) Speed of sound compared to expectations from 
soliton kinetic energy. 



In Fig. H (top), we plot the T dependence of k over 6 orders of magnitude for selected lattice sizes L with the fixed 
boundary conditions. To obtain these resuhs, we first verified that Fourier's law J — —k{T, L)VT is satisfied, by 
plotting J versus VT for various gradients. Only then, can we extract the conductivity k correctly. The dotted and 
dashed lines represent the behavior we expect from the following arguments: At high T, the quartic term for in 
the Hamiltonian H dominates over the quadratic term in q^, as can be seen from the virial theorem, for instance. 
Therefore, in this asymptotic region, the Hamiltonian has the scaling property, H /T = H' /T' under the rescalings 



[T/T'f'S'k- 



In other words, the physics at different temperatures is related by simple 



kinematical rescalings. Applying this rescaling to the thermal conductivity, we derive k = — (J)/VT — (T/T')^/^k', 
which well describes the behavior k ^ T^/* {T 3> 1) we have obtained. 

The low temperature limit can be viewed as the weak coupling limit (/? <^ 1) of H^. Near the harmonic limit, we 
expect (and verify) the speed of sound Cg and heat capacity Cy to be unity. Hence kinetic theory suggests that k ~ A, 
where A is the mean free path. The mean free path in this limit is on order of A ~ 1//3 |l2), so that k ~ 1/T. We 
also find that the size dependence of k has a power law behavior that is independent of the temperature, which is 
consistent with the previous results uqM ■ Consequently, we expect and find 



1.2L^T~ 



{T < 0.1) 
(T > 50) 



S = 0.37(3). 



(3) 



This provides a much more global picture of heat transport in the FPU model than previously known. 

The excitations in the FPU model are known to be solitonic The 'sound speed' Cg can be measured by 

examining how an initial perturbation in the center of the system propagates, as a function of lattice size L and 
temperature T. In Fig. || (bottom) wc plot this together with the velocity behavior suggested from solitons, namely, 
2T — -y/cg — 1 . The agreement is quite good aside from a weak L dependence which is not described by this 
approach. 
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FIG. 3. Behavior of VT and J as a function of T compared to expected behaviors. Solitons propagate freely below Tc ~ 



Consider now some consequences of Eq. (|l|). Kinetic theory indicates that k — XcsC^- Since Ci, ~ 1 for all 
temperatures we measured, we have A ~ n/cs- Hence, at low T, A ~ l/T, and the boundary jumps eventually 
dominate the thermal profile. The transition to this regime occurs when L ^ 2a\ — 2aK (we show below that a ^ 1), 
or using Eq. (||), Tc '-^ 2L^~^ . To delineate this behavior, we study the system at various central temperatures, T, 
with fixed relative differences AT/T = r = 0.4 in Fig. ||. For T < Tc, boundary jumps dominate Eq. (|l|), so that 
J = AT/(2a) = {r/2a)T ^ T/A. Notice here that in this regime, J should be independent of L. This is indeed 
observed in Fig. ^ (bottom). The transition also is characterized by A ^ L, which means that excitations (solitons) 
now move freely across the system without interaction. This indicates the onset of the harmonic limit which is 
characterized by a flat temperature profile. We can use Eq. (|l|) to compute how fast we approach the harmonic limit: 
VT = {r/2aK)T = Q.2T^L~^ . Hence the system approaches the harmonic limit like T^. We have explicitly measured 
this behavior for several L and confirm it. 
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A frequently used argument surmises that for fixed boundary temperatures, the behavior of k(T, L) can be obtained 
from the behavior of J by assuming that VT scales like 1/L (see for example jl^J^). Our results demonstrate that 
such an argument in general is not applicable due to the presence of boundary jumps which cannot be ignored. It 
should be emphasized that the importance of jumps is unrelated to the absence of a bulk limit 1^. 

At higher temperatures, boundary jumps compete with the heat flow inside. Here and k both behave at 
so that A is independent of T. Consequently, VT = r/(2aA + L)T = [OA/{L + 3L^)]T and J = -kVT = 

we show both VT and J as a function of T and see that the predicted behavior is 
consistent with measurements. The agreement includes not only the power law, but the prefactors as well, both for 
low and high T. 
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FIG. 4. (Top) Violation of Fourier's Law, — J/kVT, as a function of VT, showing departure from Fourier's law (dashes) at 
larger gradients. (Bottom) Temperature profiles for the three points indicated in the top figure. One can see that the profile 
has a fairly constant gradient in spite of the 20% violation of Fourier's law in (c). 



Recently we have derived a formula for the temperature profile T{x) in models where k{T) displays power law 
behavior, k ~ const. /T'^, in the temperature range of interest |,|: T{x) = Ti[l - (1 - (r2/Ti)i-^)a;/L]i/(i-T). In 
the present study, in looking at temperatures near T=8.8 in Fig. 2, k is nearly temperature independent over more 
than a decade in temperature. Hence we expect T{x) = Ti + [T2 — Ti)x/L + 0(7) to develop almost no curvature for 
very large gradients, and hence the shape of VT{x) is not a good measure of the validity of Fourier's law. In Fig. ^ 
we plot —J/kVT as a function of the measured VT. The dashed line represent the expected result when Fourier's 
law holds. One can see that as the gradient increases, noticeable systematic violations of Fourier's law develop. For 
the selected points, we plot the temperature profiles in Fig. ^ (bottom), demonstrating the nearly linear behavior. 
Similar behavior is also seen at higher temperatures, since k ~ T^/^ is a weak power law. Hence in contrast the the 
(j)'^ model, where noticeable curvature in T{x) develops when Fourier's law is broken the FPU model demonstrates 
that nearly constant temperature gradients do not imply Fourier heat flow. It should noted that one will underestimate 
the thermal conductivity if the Fourier's law is applied in this regime to obtain the conductivity. 
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FIG. 5. Quantitative study of boundary jumps 5T in our study, normalized by the mean free path A, allow the determination 
of a. The results from previous stochastic (KM) and dynamic (LLP) studies have been rescaled and are also shown. One can 
see that over an extremely large range in jumps and gradients, a simple scaling behavior exists. 
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In order to compute the coefficient a, we must analyze the behavior of the boundary jumps ET j A as a function of 
VT. In Fig. ^, we combine results from various lattice sizes and temperatures, as well as from stochastic boundary 
conditions (see for instance Figs. 1 in ref. jl^ and 2-4 in ref. Q). In Fig. ||, we plot bT jX and see that this behaves 
as VT. In addition we show the behavior extracted from stochastic boundary conditions as well as from previous 
calculations in the FPU, which are quite consistent with our results. We then extract the parameter a = 0.8(1) for 
the demons. A necessary consequence is that the temperature profiles cannot be shape invariant in the way discussed 
in 1^. While the same scaling is found for different thermostats, the coefficient a can be different. For NH demons 
(LLP) used in unh = 2.0(3) and for stochastic boundary conditions (KM) in |Q, it is much larger, as — 40(6), 

as evidenced in the large boundary jumps seen in [Q. 

We have seen that the interplay of surface and volume effects in Eq. (1) can be used to understand transport 
properties. Since the conductivity in this model has no bulk limit, the transport depends not only on the method 
used to thermalize the boundaries, but also on the boundary conditions themselves. Nevertheless, one can understand 
the temperature dependence at low and high temperatures. We have been able to map out k{T, L) and provide an 
analytical understanding in the asymptotic regimes. The FPU model also provides an example of a theory where 
Fourier's law breaks down, but the temperature profile retains a nearly constant gradient. This can be understood in 
terms of the weak temperature dependence of the conductivity. The boundary jumps also display a very systematic 
behavior over many decades in VT, independent of the types of thermostats (dynamics or stochastic) or of the lattice 
size. This dependence on VT also shows that the temperature profiles cannot be shape invariant unless the boundary 
jumps are explicitly taken into account. 
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